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1. Introduction 

An effective Hamiltonian approach is an important tool for a qualitative understand- 
ing of low energy effects in Quantum Field Theory and Statistical Mechanics [||J. The 
quantitative description involves some problems. First, it is not generally clear how to de- 
rive coupling constants of an effective Hamiltonian in terms of parameters of a microscopic 
Hamiltonian. (The latter is in common use for the study of a high energy behavior.) When 
it is not possible to relate the large scale physics to the microscopic parameters directly, 
one uses instead the coupling constants of the effective Hamiltonian as phenomenological 
parameters of the model under consideration. A related problem concerns a control over 
an accuracy of perturbative expansions. An effective Hamiltonian density is a series which 
involves an infinite set of local irrelevant fields and the corresponding perturbation the- 
ory is unrenormalizable. Thus, it is impossible to make a uniform (on the energy scale) 
estimate of an accuracy in a given order of the perturbation theory generated by the low 
energy effective Hamiltonian. 

Two dimensional exactly solvable models provide unique opportunity to refine our 
understanding of the effective Hamiltonian approach. Recently a significant progress was 
made in this field of research. In the works @ , || , || the method have been developed, 
which enables one to find an exact relation between short and long distance asymptotic 
behaviors of physically interesting quantities. As a result the coupling constants of the 
effective Hamiltonians are expressed explicitly in terms of the microscopic parameters for 
many interesting integrable models. 

This work constitutes an attempt to apply the effective Hamiltonian approach to a 
quantitative study of low-energy effects in the XXZ Heisenberg spin chain ||, 0. The 
model was chosen for several reasons. First of all, the method of the papers 0, M, M is 
suitable for calculation of the XXZ effective Hamiltonian. Next, due to more than sixty 
years of study, a huge amount of numerical data is available. Finally, the XXZ Heisenberg 
spin chain still comes to the attention of Condensed Matter physicists , || . 

In Section 2, we fix our notations and present the coupling constants of the most rele- 
vant terms in the effective XXZ Hamiltonian. In Sections 3 and 4, relying on the effective 
Hamiltonian, we develop the perturbation theory to study the low energy thermodynamics 
and finite size corrections to vacuum energies. In Section 5, we discuss the subleading long 
distance asymptotics of the two-point correlation function of the XXZ spin chain. 
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2. Low energy effective Hamiltonian 

In this work we consider the XXZ spin chain ||, @, 

T N 
Kxxz = - g E i^t+i + «i + A («+i " 1) ) , (2-1) 

(7^, <7^ and o"^ are the Pauli matrices associated with the site k of the chain. We assume 



that N is an even integer, and supplement ( |2.1| ) with the boundary conditions 

(2.2) 



„± _ e ±2iri6 ± 

°1 — e "AT+1 ' 



z z 

a l — a N+l 7 



with the real parameter < 9 < 1. Let the parameter J be positive and 

-1< A < 1 . (2.3) 



The regime ( |2.3|) is usually referred to as the "disorder regime". Notice that the unitary 
transformation, 

U+a£U=(-l) fc a£, U + ^U=(-l) fc a fc y , V + a z k U = a z k , (2.4) 

can be equivalently described by the substitution 

J^-J, A->-A. (2.5) 



In particular, the model ( [2.1|) with A = — 1 is unitary equivalent to the £TJ"(2)-invariant an- 

tiferromagnetic XXX spin chain. It will be convenient for us to parameterize the constants 

in the following manner 

A = cos(tt/3 2 ) , 

j_ 1-/3 2 _ t (2-6) 

sin(7r,9 2 ) 

with < P 2 < 1, a > 0. 

In the thermodynamic limit 

iV^oo, a^O, L = Na- fixed (iV - even) (2.7) 



the XXZ spin chain in the disorder regime ( |2.3[ ) renormalizes to the continuous quantum 
field theory - the Gaussian model M , fTO] , p] , 



lim ( H 1XZ - a x L £ ) = H Gauss . (2.1 



Here @, 

2 , „ r sinh(/3 2 t) , s 

£0 = 1-/3 2 / dt V^ — ^7T- 2.9 

iva V M ; 7o sinh(i)cosh((l-/3 2 )t) V ; 

To define the Hamiltonian of the Gaussian model let us introduce the set of operators 
satisfying the commutation relations 

[a n ,a m ] = 2n5 n+m ,o, [Q,P]=i, 

- - (2-10) 

[a n , a m ] = 2n 5 n+mfi , [Q, P] = i . 

The Heisenberg algebra admits representation in the Fock spaces, i.e. the spaces generated 
by the action of a n , ak with n, k < on the "vacuum states" \p,p) which obey the equations 



p-p = s/3 , 
p+p= (^ + m)/3 _1 , 
here m and s are arbitrary integers. It is convenient to introduce the fields 



where 



, . . _ Aixx „ v-^ 
(f>(x) = Q+ — P~iJ2 


77, 


n^O 




t(x) = Q-^P-iJ2 


a n —2g-ixn 

n 


n^Q 




;an be written in the form 




n G auss = J Q ^{ T + T}, 


T(t) — • (8 rh) 2 ■ 


n 2 


^ ' ~ 4 ' ^° x ^> ■ 


QL 2 ' 


T('r) — • (f) rh\ 2 ■ 


7T 2 


\ ' ~ 4 ' \ x ™> ' 


6L 2 ' 
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(2-11) 



an\p,P) = a n \p,p) = for n>0, 
P|p,p) =p\p,p), P\p,P) =P\p,p) ■ 
Then 

Hcauss = J;{ P2+P2 ~V2 + 2^ {a-nttn + a- n a n ) } . (2.12) 

n>0 

This operator acts in the Hilbert space 

it = @p,p~Fp,p 5 (2-13) 

with 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



To describe the bosonization of the matrix er?, we define the fields 



(p(x) =<j>-<t> 

ip(x) = <fi + (ft 



(2.18) 



In the leading order in a [pja, 



E- 
a k — a 2 



°k ~ 



2^W 



(2.19) 



2tv(3 



d x (p(x) 



with x = ka. To give a precise meaning to the constant F one needs to specify the normal- 
ization of the exponential fields. We will normalize exponential operators in accordance 
with the short distance expansion 



D icnp 



(x) e 



- ia *{x' 



e %0 " f, (x)e~ t0 " f> (x') 






CXO lAJ tij 







(2.20) 



for an arbitrary a. In this normalization the constant F = F((3 2 ) was found in Ref. [13 

1 



2(l-,3 2 ) 2 

exp 



W P 2 \ 10 

1 \2-2B^J 



x dt ( sinh(/3 2 t) 



X 



o t V sinh(t) cosh((l - (3 2 )t) 



(3 2 e~ 2t ) } . 



(2.21) 



The Hamiltonians ( |2.1| ) and ( |2.1(j| ) coincide in the lowest nontrivial order of the lattice 
parameter a (|2.8|). In the next orders irrelevant fields contribute to the density of Hxxz- 
These irrelevant fields should be local with respect to the spin fields ( |2.19| ) and commute 
with the operator 

1 N 



(2.22) 



s=l 



A simple analysis shows that the most important corrections have the form |TT 

Hxxz=a- 1 L^ + | L ^{T + T-a 2 (A + TT + A_(T 2 + T 2 )) + 

a^" 2 A cos (<p/0) +...}. 



(2.23) 



In the case of a%, Eq.(2.19) gives the leading order in a for < [3 2 < 



2 " I 

2" 



Here, the field T 2 (T 2 ) is defined as a regular part of the operator product expansion of two 
operators T (T). The dots in ( |2.23|) and below mean a contribution of higher dimensional 
local counterterms. In the Appendix we describe the calculation of the coupling constants 
in the low energy effective Hamiltonian ( |2.23| ) , 



A 



4T(/r 



T(l-(3 



r (i + idbO 



i 



A+ = 27 tan (2 



TV 



A_ 







L2v^r(i + ^ 
)> 

*-\2-2B 2 ) F V2-2/3 2 ) 



) 



-2(3 2 

yr 3 ( 



(2.24) 



i2-r(^i,)r3(^) 



For the XXX spin chain (A = — 1) 



A_ 



0, 

4tt 



(2.25) 



We discuss how to treat "A-term" from ( |2.23| ) at the limit /3 2 — ► 1 in the next section. 

It is important to note that except for the "A_-term" , the Hamiltonian ( |2.23|) coincides 
with low energy effective Hamiltonian for the sausage model with the value of the #-angle 
equals to tx || .□ The fields T 2 and T 2 have Lorenz spins 4 and —4 respectively. Therefore, 
the "A_-term" destroys the rotational symmetry of the Gaussian model in the infinite 



volume and appears because of the lattice nature of the Hamiltonian (|2.1|) . 



3. Temperature and magnetic field corrections to specific free energy 



The effective Hamiltonian ( |2.23|) can be applied to the study of infrared corrections 
of physically interesting quantities in the XXZ spin chain. Here we deal with the low 
temperature thermodynamics. Let us consider the model (|2.1| ) at the low temperature 
T (T/J < 1 ) and in the weak magnetic field h ( \h/J\ < 1). With the result (gjgjj ) one 



V.A. Fateev et al. use the notation A and M which are related with the parameters (3 2 and 
a as following, 



A = 1 - /r 



M = Aa 



can develop the perturbation expansion in the small parameter a of the specific free energy 
fxxz, 



Notice that "A±-terms" in ( |2.23| ) contribute to the first order of the perturbation series 
and lead to a 3 -corrections. The term with A does not have diagonal matrix elements 
and its contribution appears only in the second order. This implies a correction which is 
proportional too? . Therefore the leading corrections to the specific free energy differ 
significantly in the domains < (3 2 < | and | < (3 2 < 1 and have the form 

fxxz=£o + af Gauss -^ (a+((T» 2 + 2A_((T 2 )))+... for < /3 2 < ? , 

fxxz = £o + a fa 
X 2 J 



rauss 
-■■>> rL 



16tt 2 



dx / dr (( e^ /f3 (x, r) e -^ //3 (0, 0) )) + 
(I Jo 



for - < P 2 < 1 . 



(3.2) 



Here we use the notation 

p ^ JGauss —— 'TV- 



n 



p T Gauss \ 2(3T V / 



<(0» = Tr w 



„— y HGauss+2pr (P~ P) O 



Tr 



n 



- TjT Hqo 



'+ 2/3T ( " ") 



and 



<p(x, r) = 0(x + it) + 4>{x — it) 



The calculation of ( |3.2|) is straightforward. We will not discuss here the most general 
expressions but restrict our attention to the case T 3> L~ l . Under this condition the 
system develops an effective correlation length and all finite size effects are suppressed by 
the exponentially small factor e~ nLT . Omitting the finite size terms, one obtains 

/x„ = ft + 4 TO3 -/ 1 ( g ^)-(2™)'T«{A + i?(* J+2A. /,(*)} + .... (JU.) 



for < (3 2 < |. If | < (3 2 < 1, the leading correction comes from the second order 
perturbation theory term 



fxxz =S Q + ^aT 2 I l( ^)- 



4aT 2 ^sin(^) 



aT^T(l + ^)^--' 



6 



H 



ih 



8tt(3T 



(3.4) 



+ 



The auxiliary functions in ( |3.3| 



read explicitly 



h( P ) = 
h{p) - 
Hi(p) 



P "24 



P 



P 

T 



4/3 4 - 17/3 2 + 4 



960/3 2 

r(/r 2 )r(i-2/r 2 ) r(r 2 + 2p/3- 1 ) 



(3.5) 



/3 4 r(l-/3" 2 ) r(l-/3" 2 + 2p/3- 1 ) 

and A± = A±(/3 2 ) are given by (|2.24| ). Notice that hip), Ia(p) coincide with the vacuum 
eigenvalues of the first local integrals of motion in the quantum KdV theory |14[] . At the 
same time H\ (p) is a vacuum eigenvalue of the first "dual unlocal" integral of motion |TjJ . 
If (3 2 = x 2 2 n (n = 1,2, ...), the first order contribution of the "A+-term" interferes 
with the 2n order of the perturbation expansion generated by the "A-term" . The poles of 
the coupling constant A+ ( |2.24| ) at these points indicate this phenomena. The interference 

I the specific free energy 



produces a logarithmic contribution. For example, at /3 2 
behaves as 

2 
3^3 a 96v 

11 . / x „ / , ,1 3z/i 



JXXZ \q2 



±- (l6n 2 T 2 + 9h 2 ) (l + |a 2 T 2 + -^ (16tt 2 T 2 + 9h 2 ) 
ii-l 0g ( a r/12)- S e(^(I^)) 



1927T 2 



X 



(3.6) 
where ip(t) = d t log (r(t)). 

If (3 2 ~ 1 , Eq. ( |3.4|) defines the leading asymptotic behavior in the very narrow domain 
of temperature and magnetic field, 

log (^' log (^» >> T 1 



P 2 

In the limit /3 2 — * 1 the domain of validity of ( |3.4| ) disappears completely. Hence the case 
of the XXX spin chain requires an additional analysis similar to the one of Refs. @ , . 
Let us consider the Gaussian model with f3 2 = 1. The fields 



Jo 
Jo 



\d x c 
\d x c 



J± 

J± 



,±io 



,=F*9 



generate the right and left current algebras at the level k 
coincides with the Wess-Zumino-Witten (WZW) model, 



(3.7) 



1, and the Gaussian model 



H 



Gc 



p 



H 



WZW 



(3.8) 



In the vicinity of the point /3 2 = 1 it is convenient to rewrite the effective XXZ Hamiltonian 
as the marginal current- current perturbation of the WZW Hamiltonian [[ll| , 

Hxxz = a~ l L£ Q + H wzw + J ^ {g\\ JqJq + y (J+J- + J-J+) +■■■}■ (3.9) 

Here g» > and \g±\ < g\\ are small running coupling constants. The corresponding 
Renormalization Group (RG) equations are known exactly [Q, |l7[j . Under a suitable 
diffeomorphism of gn and g±_ (i.e. the renormalization scheme choice), one can set up 



dr 
dg± 



2-<7|| 
2 0n0± 



(3.10) 



dr 2 - 5r|| 

Here r is the RG length scale. These equations are solved as 



<7H = 2(l-/3 



2, 1 + g 
1-9 



^ = -4(l-/3 2 ) *' 



1-9 



with 



r 



4 "^ 



(3.11) 
(3.12) 



while ro is a r-independent constant. 

The correction to the specific free energy, descended from the "A-term" of ( 2.23|) , 
admits a power series expansion, 

anT 2 

- i 1 + «io 0|| + "20 PiT + «02 01 + U[g") \- 

(3.13) 



/xxz = £o — [l + wio 0|| + w 20 g\ + u 02 g\ + 0(# 3 )}- 

|l + u 10 fir|| + U20 tf 2 + ^02 #1 + 0(5T 3 )| + ... . 



ah 2 r 



To find values of the first coefficients u^j , Vkj we should compare the expansion of ( |3.13j ) 
and (|3.4| ) in power series of 1 — (3 2 for a fixed value of the parameter q. According to ( |3.12|) , 
the latter can be chosen as a solution of the algebraic equation 



q (1 - q)~P 



aTe-i^T(l + ^) 



r (!+2^) 
A simple calculation leads to the result 

ixaT 2 r 3 



r(r 2 + 4^r) 



r(2-/3- 2 + ^w) 



(3.14) 



fxxz = £o-^-{i + l 9l] gi + 0( 9 *)}- 



6 

lftr l 2 4 4 8 32 vy ; 



(3.15) 



+ ... 



Now we can take the limit 1 — > 1, 

#u ->• g , g±->-g , (3.16) 



of Eqs. (|3.11|) , ( |3.14f) , (|3.15|) and obtain the low temperature behavior of the specific free 



energy of the XXX spin chain, 

/jc „=-2/ I o g (2)-g{l + | 9 « + to«)}- :[ ^ 7 {l + | + | 9 « + te«)}- 

y/3 f ^ 4 h 2 T 2 37r 2 T 4 ^ , ! 
16^1 64^ + ~ + -5- +O0/) } + •••• 

(3.17) 

The function g = g(T, h) in ( |3.17|) solves the equation 

< 7 - 1 + ilog(<7) = -Ke(^(l + ^))+log(v / 2^e^J/T) , (3.18) 



with ?/>(£) = <9 t log (r(t)). In the formulas ( |3.17| ), ( |3.18j ) we use the original parameter J 



of the Hamiltonian (|2.1| ) with A = 1. We also include the correction arising from the 
"A_-term" ( |3.3| ). From ( |3.18| ) one has the formula for the magnetic susceptibility, 

x = -4d 2 h f xxx \ h=Q , (3.19) 

as a function of temperature, 

Now the running constant depends on temperature only, g = g(T), 

g- 1 + l - log(<7) = log (V^e^J/T) 

and 7 = 0.577216... is the Euler constant. The leading T-dependence (the first order in g) 
of x(T) was obtained previously in jn|. In Fig. 1, ( |3.20| ) is compared against the result of 



numerical solution of the Thermodynamic Bethe Ansatz equations (T^ 




T/2J 



Fig. 1. Magnetic susceptibility of the XXX antiferromagnetic spin chain. The 
continuous line represents ( |3.20| ). The bullets were obtained by the numerical 
solution of the Thermodynamic Bethe Ansatz equations [|18] . 



4. Finite size corrections to vacuum energies 

In this Section we consider the XXZ spin chain at zero magnetic field and temperature. 
If T = 0, the function fxxz ( ft-l| ) coincides with the specific ground state energy E xxz /N 
and Eqs.(|3.2|) are useful to study finite size effects. In fact, the calculation can be easily 
generalized to obtain the finite size corrections to an arbitrary vacuum eigenvalue (not 
only ground state) E xxz in the sector with a given value s of the total spin S z ( |2.22|) . 
Let us parameterize E xxz /N by the form, 



E 



XXZ 



/N = £ + 



2n 

^N 2 



6(s,6,N) 



(4.1) 



Then, 



t s 2 f3 2 e 2 

5 = — £— + 



1 



+ S cos + S T + ... . 



(4.2) 



2 ' 2(3 2 12 
The function St gives the first order contribution of the "A±-terms" in ( |2.23|) . It reads 

10 



explicitly^ 



rT 47T 2 f A T ,s(3 9 , T f sj3 9 s , / T ,s3 6, T ,sj3 9 N \ i , An . 

The term S cos in ( |4.2|) arises from the second order perturbation theory generated by the 
field cos(<p//3), 



2 n /27Tn 



0FT(1+; 


P 2 ^ 


AT(l + 7 


l 

-2/3 2 



7^~ 



*(¥+£>*(¥-£>■ ^ 



The auxiliary functions in (fO|) , (£Oj) are given by ( |3.5| ) . For /3 2 sufficiently close to unity 
the finite size corrections require an additional RG consideration. In a similar way to 
( P-15|) , an appropriate RG analysis leads to the following expression for (| 



8(s,9,N) = 8 uu + 8 x + ... . 
The function 5 RG is a RG improved contribution of the "A-term" . 



(4.5) 



;RG 



2124 48 32i vy; ' 

where the running constants are given by (|3.11] ) with 

v^e-ir(i + ^) 



(4.6) 



q (1-q) 



The function 



r(/3" 2 + s + 9[5~ 2 ) F((3- 2 + s - 9(3~ 2 ) 
r(2 - P~ 2 + s + 9[5~ 2 ) r(2 - (3~ 2 + s - 9(3~ 2 ) 



5{s,9,N)-5{0,0,N) 



(4.7) 



was calculated in |Tl]] by the numerical solution of the Algebraic Bethe Ansatz equations 



for f3 2 = | and some vacuum energies. In Table 1 the numerical results are compared 
against the derived formula. The following comment is appropriate here. The function 



In writing ( |4.2| ), we assume that j3 2 ^ 2 n+i > n = 2,3,... (see comment to the equation 

y)). 



ii 



p 



6 i s 

2/3 + 2/3 



as well as the right hand side of (|4.7|) , has the pole at s = 0, 6 = 1. The reason of 
the singularity is rather simple; We have developed the perturbation theory starting with 
the Fock vacuum states \p,p) Q2.ll] ), with 

6 s 
' P ~ 2/3 ~ 2/3 ' 

But for s = and 6 sufficiently close to unity, the contribution of the vacuum state with 
p = p = -=3- to the true ground state of the XXZ spin chain becomes significant. It seems 
reasonable that the disagreement between the numerical data and obtained asymptotics 
for the case s = 0, 6 = h in Table 1 is due to this effect. 

The finite size corrections for the XXX spin chain can be derived by taking the limit 
(ETIBD of Eqs.(]4l|, OOP . In particular, 



7T 



v/3 



2/)2 



;RG\ 



l/3 2 =l 



8N 2 

1 
~12 



{ S 4 + ^ 4 + 6 S 2 ^ 



l+o^ + 



s z - 6 Z + 



9_,9_ 
2 4 



20 

2 7 

g 3 



32 



16 l y 



r7 3 }+ (4-8) 



?{ 1 + i + f} + <^) 



with 



The leading terms in ( |4.8| ) for ^ 
IT 



g- 1 + i log(^) = log (2K-h*N) -^{l + s + 6)+iP(l + s-6)y 

0, 1 are in agreement with the results of the work 



0, s 





s = 0, 6 = \ 


s = l, = 


s = l,0=i 


N 


BA 


RG 


BA 


RG 


BA 


RG 


8 


0.15519 


0.15467 


0.39793 


0.395... 


0.43606 


0.433... 


16 


0.15217 


0.15185 


0.40594 


0.40555 


0.44460 


0.44406 


32 


0.15104 


0.15081 


0.41063 


0.41060 


0.44897 


0.44883 


64 


0.15055 


0.15037 


0.41327 


0.41327 


0.45128 


0.45121 


128 


0.15030 


0.15015 


0.41474 


0.41474 


0.45254 


0.45250 


256 


0.15017 


0.15003 


0.41557 


0.41557 


0.45324 


0.45321 


oo 


0.15000 


0.41600 


0.45416 



Table 1. The function 5(s, 6, N) - 5(0, 0, N) for f3 2 = |. The "BA" columns 
were obtained by means of the numerical solution of the Algebraic Bethe Ansatz 
equations |I[. The "RG" columns follow from (p|), flOj), (fO|). 
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5. Subleading asymptotics of correlation functions 

As well as the Hamiltonian density (|2.23| ) , the formulas ( |2.19|) give only leading terms 
of the expansions in local quantum fields of the spin operators. The leading terms allow 
one to obtain the leading asymptotic behavior of correlation functions of the spin opera- 
tors. Subleading asymptotics are the result of the two effects: next terms in the expansion 
( P-19|) and XXZ corrections of the Gaussian ground state, which defined by the effective 
Hamiltonian ( |2.23| ). Therefore, to derive the subleading asymptotics of the correlation 
functions systematically, one needs to know next terms in the expansion ( |2.19| ). Unfor- 
tunately, this problem admits only a qualitative analysis at the moment. Nevertheless, 
some interesting quantitative predictions for the subleading asymptotics can be obtained 
based on the effective Hamiltonian (|2.23j ) only. As an example we consider the equal time 
correlator 

W(n) = : r , n = l,2... , (5.1) 

{vac | vac) 



for the infinite spin chain (L = oo). Simple arguments |16[] show that the most important 



correction of the bosonization formula for a k , n has the form 

4 +n = a^ sj^e ±l ^(x) (l±i(-l) n aAd x <p(x) + -) {x = na) . (5.2) 

The exact value of the coefficient A is not known nowadays. It is instructive to note that 
due to the sign factor (— l) n in ( |5.2| ) a linear contribution of the "A-term" in W(n) is 
canceled out exactly for odd integer n and the subleading asymptotic is determined by the 
effective Hamiltonian only. We conclude that for | < (3 2 < 1, 

W(n) =Fn- p2 {\-Cn~i 1 - 2/3A ((-l) n + l) n" 1 + ... ] , as n -»• +oo . (5.3) 

Here the constant F is given by ( |2.21|) , A = A(j3) is not known, while the constant 
C = C((3 2 ) can be obtained by the second order perturbation theory based on the effective 
Hamiltonian ( |2.23| ). The calculation leads to the formula 



c r2</3 



r 2 (i-/3- 2 ) 



r r (! + 2^) 


- 


[2^T(l + ^ 


oJ 



^" 4 / 2^ 



sin 



(27T/3- 2 ) 



(5.4) 



P i 1 1 n-2\ ilt^ 



(1_ /?2) ( 2-/32) 
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#(r*)-n\-r 2 )) 



with ip f (t) = df log (r(t)). Eq. (|5.3|) defines the subleading asymptotic behavior of the 
correlation function for 

In order to study the limit /3 2 — > 1, it is useful to rearrange the expansion of W(n). To do 
it, let us introduce g\\ and g± by Eqs. (|3.11| ), while 



q (l-q)' 



e^rfi + A) 



and 7 is the Euler constant. Using (|5~3|), (|5.4|), one obtains 



7P"" 



(5.5) 



^ (n) = M m {i -£ 



<7i 4C(3) + 3 3 , 164C(3)-67 2 



8 16 



192 



»ll 



+ 



384 



^4 

2n 



((-iy + i + (g))+0(n- 2 )} 



9\\9± + 0(g*)- 



(5.6) 



where .A = 4/5A„ 2 _ ] and £(s) is the Riemann zeta function. The non-perturbative func- 
tion Z in (5.6) reads 



Z(g\\,g 



jl, 



9±\ 



L2e2 ^|| 

We can derive the asymptotic of the correlation function of the XXX spin chain by taking 
the limit g\\ —> g , g± —* —g • Finally one haaa 



W(n) 



l/3 2 =l 




7T 3 n^fg 



i r. 3 2 , 156C(3)-73 3 , _, 4 , 



1 g* + 

16 y 



384 



^ + 0(g* 



A 

In 



((-iy + l + 0(g))+0(n- 2 )} 



and 



g- 1 + ~ log(</) = log (2v / 2^e^+ 1 n) . 



(5.7) 



The correlation function VF(n)L 2 _ ] (1 < w < 30) was calculated numerically in the 
works [ZT|, ||22| B In Table 2 the numerical data are compared against ( |5.7| ) for odd integers 



The coefficient y -\ in Eq.(5.7) was derived previously in [20 based on the explicit form of 
the constant F ( |2.2l| ). I am grateful to I. Affleck for communications concerning this calculation. 



K.A. Hallberg et al. [22 used the notations u{n) = W{n)/A 



/3 2 = l' 



14 



0.16 



0.14 



0.12 



n W(n)/4 




o.os 



Fig. 2. The correlation function j 



{vac\a*a* n \vac) 



for the XXX spin chain. 



4 ( vac | vac) 

The bullets (see Table 2) were obtained in [E21. The continuous line I follows 



from (|5.7| ) for odd integers n. The line II represents ( p.7| ) for even n ( A. = 0.4 ). 



n. Unfortunately, we can not directly compare the numerics for even integers n, since the 
value of the constant A in (|5.7|) is not known. Fitting of the data gives 



-4^0.4 



(5. 



The corresponding plots are presented on Fig. 2. 

Finally we note that there are two different approaches to exact calculation of cor- 



relation functions like ( |5.1| ) at the moment. In the first one the correlation functions are 
expressed in terms of Fredholm determinants @, [23], [24]. Recently M. Jimbo and T. 
Miwa |2!| found another representation for the correlation functions in terms of n-fold 
integrals. It would be interesting to derive the large n asymptotics (|5.3| ), ( |5.7|) via these 
approaches. 
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n 


NUM 


RG 


1 


0.1477 


0.1239 


2 


0.1214 




3 


0.1510 


0.1427 


4 


0.1384 




5 


0.1541 


0.1505 


6 


0.1463 




7 


0.1567 


0.1554 


8 


0.1513 




9 


0.1596 


0.1589 


10 


0.1550 




11 


0.1620 


0.1616 


12 


0.1581 




13 


0.1641 


0.1639 


14 


0.1606 




15 


0.1659 


0.1658 


16 


0.1628 




17 


0.1676 


0.1674 


18 


0.1646 




19 


0.1689 


0.1689 


20 


0.1661 




21 


0.1700 


0.1702 


22 


0.1674 




23 


0.1712 


0.1713 


24 


0.1687 




25 


0.1723 


0.1724 


26 


0.1699 




27 


0.1734 


0.1733 


28 


0.1710 




29 


0.1746 


0.1742 


30 


0.1722 





Table 2. The correlation function j 



'k + n 



chain. The column "NUM" were obtained in |22 
from ( |5.1| ) and ( |5.7|) for odd integers n. 



i for the XXX spin 

( vac | vac ) r 

The column "RG" follows 
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6. Appendix 

The calculation of the coupling constants in the effective Hamiltonian (|2.23| ) is based 
on the technique developed in the works 0, 0, ||. We discuss it very briefly. 

Let us consider the XXZ spin chain ( |2.1| ) with infinite number of sites N in the 
magnetic field h ( |3.1|) . The specific ground state energy 



reads [[T2| 



8(h) 



h 



lim E 



(0) 
XXZ 



(h)/N 



e{h) = -- + j 



A 



da 



sin(7r/3^ 



e a 



_ A 2n cosh(2ct) + cos(7r/3 2 
where e(a) satisfies the Yang- Yang equation 

" A da' sin(27r/5 2 ) e(a') 

_ A tv cosh(2a — 2a') — cos(2irf3 2 ) 



(6.1) 



(6.2) 



e (a) 



h 



4sin 2 (7r/3 2 



J cosh(ct) + cos(2n/3 2 ) 



(6.3) 



Solving 
£{h) = So 



e(±A) = . 
by the Wiener-Hopf method, we obtain for the weak magnetic field 

"leW (i + ^l^l^" 4 + (^+ + ^-)(^) 2 + -) (IV^I«i), (6-4) 



where the constant a and £q are given by (|2.6|), (pT9|), and 



K 



«+ 



«. 



r 2 (/?- 2 )tan(7r/3- 2 ) 
2/3 2 r 2 (i+/3" 2 ) 

1 , TX x 

= 64^ ^ ( 2^2^ 



r (2-2S 2 ) 



V^ r (2^W) 



(6.5) 



1 r (2-2/32) r (2-2/32) 



,2-2/3 
~3/32 



i92tt w 3/j , )r 3 ( 1 ., ') 

1 l 2-2/32^ J- { 2 -2/3 2 ) 



Now, we should compare (|6.4| ) with the result of the perturbation theory based on the 
effective Hamiltonian (|2.23|) . This allows to extract the values of the constant A and the 
combination A+ + 2A_. The "A+" and "A_-terms" in the effective Hamiltonian ( |2.23| ) differ 
essentially. The operator TT has Lorenz spin 0, while T 2 (T 2 ) has Lorenz spin 4 (—4). 
Hence, the A+-term preserves the Lorenz invariance, whereas the "A_-term" destroys the 
rotational symmetry. A careful analysis of the expansion Q6.4J ) based on this simple obser- 
vation makes it possible to split contributions of the operators TT and T 2 + T 2 and leads 
to Q2T241) . 
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